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■ ABSTRACT 

\^ \ Recent results indicate a correlation between nuclear radio-loudness of active galaxies 

(-H . and their central stellar surface-brightness profiles, in that 'core' galaxies (with inner 

^ ^ logarithmic slope 7 ^ 0.3) are significantly more radio loud than 'power-law' galaxies 

I I (7 ^ 0.5). This connection, which indicates possible links between radio-loudness 

2 ' and galaxy formation history (e.g. through black hole spin) has so far only been 

" ■ ' confirmed for a radio-selected sample of galaxies. Furthermore, it has since been shown 

that the Nuker law, which was used to parameterise the brightness profiles in these 
studies, gives a poor description of the brightness profile, with its parameters varying 
systematically with the radial fitted extent of the profile. Here, we present an analysis 
of the central surface brightness profiles of the active galaxies of Hubble Type T ^ 3, 
that were identified by the optically-selected Palomar spectroscopic survey of nearby 
galaxies. We fit the brightness profiles using Sersic, Core-Sersic and, where necessary, 
Double-Sersic models, which we fit to the semi-major axis brightness profiles extracted 
from high resolution images of the galaxies from the Hubble Space Telescope (HST). 
We use these fits to classify the galaxies as 'Core', 'Sersic' or 'Double-Sersic'. We 
">) ' compare the properties of the Active Galactic Nuclei (AGNs) and their host galaxies 

with this classification, and we recover the already established trend for Core galaxies 
to be more luminous and contain a higher-mass supermassive black hole. Defining the 
radio-loudness of an AGN as the ratio of the nuclear radio luminosity to [Oiii] line 
luminosity, which allows us to include most of the AGN in our sample and prevents a 
bias against dim nuclei that are harder to extract from the brightness profiles, we find 
that AGN hosted in Core galaxies are generally more radio-loud than those hosted 
, in Sersic galaxies, although there is a large overlap between the two subsamples. The 

correlation between radio-loudness and brightness profile can partly be explained by 
a correlation between radio-loudness and black hole mass. Additionally, there is a 
significant (99 per cent confidence) partial correlation between radio-loudness and 
the Core/Sersic classification of the host galaxy, which lends support to the previous 
results based on the radio-selected sample, although it is possible that this partial 
correlation arises because AGN in core galaxies tend to have a lower accretion rate as 
well as a higher central black hole mass. 
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1 INTRODUCTION 

Before the Hubble Space Telescope {HST) was launched, ob- 
servations of the central regions of galaxies were limited by 
atmospheric seeing, which prevented resolutions better than 
^ 1" from being achieved. This made it difficult to deter- 
mine how the brightness profiles of these galaxies behaved at 
small radii. However, high resolution images taken with the 



HST (Cra ne et al.lll993 : iFerrarese et al. I I1994I : iForbes et al.l 
11995) demonstrated that the surface brightness profiles of el- 
liptical galaxies and disk galaxy bulges continue to increase 
up to the resolution limit, confirming the re sults of ground- 
based observations (jKormendvll 19851 : iLauerl 1985) that con- 
cluded that the profiles of most eUiptical galaxies cannot 
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be described by models with flat, isothermal cores such as 
described bv lKind l| 19661 ). 

Early studies using HST suggested that the central 
brightness profiles can be divided into two distinct classes — 
in "power-law" galajcies, the logarithmic slope of the bright- 
ness profile remained steep towards the resolution limit of 
HST, while in "core" galaxies the slope flattened signifi- 
cantl y. To describe these brightness profiles in the central re- 
gions, [Lauer|etlal] (|l995i ) introduced a form of broken power 
law, the so-called "Nuker-law", with the logarithmic slope 
in the inner region, 7, being flatter than in the outer region. 
These classes then correspond to a dichotomy in the val- 
ues of 7, with power-law galajcies hav i ng 7 ^ 0.5 and core 
galaxies having 7 ^ 0.3. iFaber et al.l (|l997l ) demonstrated 
that there were no galaxies in t he range 0.3 < 7 < 0.5, al - 
though subsequent studies (e.g. iRavindranat h et al.l I2OOII ) 
have found a small number of such "intermediate" galaxies. 

This dichotomy in brightness profiles has been linked 
to global galaxy properties, for example core galaxies tend 
to be very luminous, slowly rotating galaxies with boxy 
isophotes, while power-law galaxies tend to be les s luminous, 
rapid l y rotating galaxie s with disky isophotes (Ijaffe et al.l 
I1994I : iFaber et al.l Il997l ). For galaxies that contain active 
galactic nuclei (AGNs) some of these properties of the host 
galaxy also correlate with properties of the AGN, for ex- 
ample the mass of the central supermassive black hole 
(SMBH) correlat es strongly with the velocity dispersion of 
the h ost galaxy l|Ferrarese fc Merritt]|2000l : iGebhardt et all 
I2OOOD . suggesting that the evolution of the host galaxy and 
AGN are interdependent. Investigating how the brightness 
profiles of the central regions of these host galaxies are re- 
lated to the properties of the AGN is important in under- 
standing how they influence one another and how t hey form 
and evolve together, for example lFaber et al.l (|l997l ) propose 
a model whereby a merger between two galaxies can pro- 
duce a core galaxy as their black holes sink to the centre, 
ejecti ng stars in the centra l region and scouring out a core, 
while Ivan der Marell (Il999l ) suggests that the cusp slopes of 
core and power-law galaxies can be reproduced by adiabatic 
growth of the black hole in an isothermal core. 

More recent s tudies (e.g. iGraham et al.l l2003l : 
iFerrarese et~al] l2006a^ have shown that the Nuker model 
does not give an accurate parameterization of galaxy 
surface brightness proflles. It was only intended to model 
the central regions and so cannot describe the brightness 
profiles at larger radii, where the outer power law of the 
Nuker model fails to properly fit the cu rvature of the 
bright ness profiles in most galaxies, and iGraham et al.l 
l|2003l) demonstrate that the parameters of the Nuker model 
are sensitive to the radial extent of the profile to which it 
is fitted. However, the brightness profiles of "power-law" 
galaxies are well fitted by either a Sersic law over their 
entire radial extent, or by a Sersic law modified with the 
addition of a second Sersic component, the latter meant to 
represent an inner compact stellar nucleus (generally below 
a few tenths of an arcsec). Meanwhile, "core" galaxies 
show a central defici t with respect to thei r outer Sersic 
profile l|Truiillo et al.l |2004| : IFerrarese et al. I |2006b). This 
provides an alternative classification to the core/power 
law dichotomy based on inner logarithmic slope, in which 
galaxies can be separated into those that are well fitted by 
a Sersic law at all radii, those that require an additional 



inner Sersic component to account for an inner compact 
stellar nucleus and those that show a central deficit from 
their outer Sersic profile. The profiles of galaxies that show 
a central deficit can be parameterized by a "Core-Sersic" 
model (introduced by IGraham et alll2003l ) that combines a 
Sersic profile at large radii with a power law in the central 
regions (see e.g. lGraham|[20o3 : ICote et al.|[2007l ). 

Recently, C apetti fc Balmaverde publis hed a series 
of three papers llCapetti fc Balmaverdel |2005| & 2006 and 
iBalmaverde fc Capetti 20061 ) in which thev investigate the 
core/power-law dichotomy, based on the Nuker law, in two 
radio-selected samples of galaxies that are likely to host an 
AGN, requiring a minimum radio fiux of ~ 1 mjy. Us- 
ing new Nuker fits to the data as well as fits from the 
literature, they showed that radio-loud AGNs (defined us- 
ing the ratio of 5 GHz to B-band luminosity densities) are 
hosted by core galaxies, while radio-quiet AGN s are h osted 
by power-law galaxies. ICapetti fc Balmaverdel (|2007l ) also 
showed that early-type galaxies hosting Seyfert nuclei show 
power-law type central brightness profiles. This correlation 
of AGN radio-loudness with brightness profile offers impor- 
tant clues to the origin of the radio-quiet/radio- loud di- 
chotomy in AGN, as well as the link between the central 
black hole and the host galaxy. For example, if the merger 
history governs whether a galaxy is core or power-law, we 
might speculate that core galaxies host rapidly spinning 
black holes resulting from a major merger, and that radio- 
loudness is linked to spin. Alternatively, the different stel- 
lar density profiles in core and power-law galaxies may be 
linked to different environments which may enhance or in- 
hibit jet production, e.g. lower or higher gas densities respec- 
tively, and /or d ifferent types of hot or cold accretion (e.g. 
iBest et al.ll2005h . Regardless of its physical interpretation, 
the radio- loudness/brightness-proflle connection may shed 
useful light on the connection between radio-loud AGN and 
giant elliptical hosts, since core type galaxies are almost all 
giant ellipticals. 

The implications of the radio- loudness/brightness- 
profile connection in AGN make it worthy of further study. 
In particular, it is important to identify whether the link 
persists in an optically selected sample of AGN. This is be- 
cause the Capetti fc Balmaverde sample was radio-selected, 
so that it is still possible that there exists a population of 
AGN which inhabit core galaxies but are sufficiently radio- 
quiet that they were not selected. Clearly it is technically 
very difficult to measure central brightness proflles for op- 
tically luminous AGN, but in any case the Capetti fc Bal- 
maverde sample is dominated by AGN with relatively low 
optical luminosities, i.e. predominantly optically classified 
as Low lonisation Nuclear Emission Regions (LINERs) or 
other Low Luminosity Active Galactic Nuclei (LLAGNs). 

In o rder to search for such AGN in nearby galaxies, 
et al. I il995) conducted an optical spectroscopic survey 
of the nuclei of bright, nearby galaxies in the Northern hemi- 
sphere, the so-called Palomar spectroscopic survey of nearby 
galaxies. Their survey included 197 galaxies that are con- 
firmed to contain an AGN, as measured by their emission 
line ratios. This gives us an optically selected sample of low- 
luminosity AGNs, in contrast to Capetti fc Balmaverde's 
radio-selected sample. In this paper, we use available HST 
images to analyse the brightness proflles of the Palomar sur- 
vey active galaxies, using the Sersic/Core-Sersic classiflca- 
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tion scheme rather than the core/power-law classification 
obtained from the Nuker model. We then use these results 
to investigate how the radio properties of these AGNs are 
related to their Sersic/Core-Sersic classification, which al- 
lows us to show whether Capetti & Balmaverde's results 
are reproduced in a relatively complete sample of optically- 
selected AGN using this alternative classification scheme. 



2 SAMPLE AND DATA REDUCTION 

2.1 The Palomar spectroscopic survey of nearby 
galaxy nuclei 

To reliably detect the presence of a galaxy core it must 
be nearby so that a core could easily be resolved by the 
HST. The Palomar sp ectroscopic survey of nearby galaxy 
nuclei JHo et al.1 [l995l l uses a fiux-limited sample of 488 
bright Northern galaxies, with Bt ^ 12.5 mag and decli- 
nations (5 > 0°, and naturally satisfies our requirement for 
nearby galaxies. Spectroscopic obse rvations of these ga laxies 
show that they include 197 AGNs. iNaear et aP l|2005l ) con- 
ducted a high resolution radio survey of all of these AGNs, 
which will enable us to compare their radio properties to 
their Core/Sersic classification. For spiral galaxies it is the 
brightness profile of the bulge component that determines 
this classification, however late-type spiral galaxies have rel- 
atively small bulges compared to the disk component, which 
will make it more difficult to cl assify them. Therefo re we take 
as our sample the 150 AGNs in lNagar et al.l l|2005l ) for which 
the host galaxy is of Hubble type T ^ 3. 

IHo et all |l993) list various parameters for the galax- 
ies in the Palomar sample, such as distances, Hubble 
types, AGN classification and bulge and total B-band ab- 
solute magnitudes, which we use in this paper. However, 
where galaxies hav e mo re recent distance estimates listed in 
iTonrv et al] (l200lj) and iPatureF et al. (200^), we have used 
these values, and updated the absolute magnitudes and lu- 
minosities accordingly. The [Om] and H/3 narrow line lumi- 
nosities for the AGN were obtained using the combination 
of u nrcddcned line ratio data and Ha fiuxes in Table 2 of 
I Ho e t al. (1993)- We d ereddened these fluxes by applyin g 
the extinction curve of ICardelH. Clayton fc M^thii (l989f) . 
assuming the Balmer decrement to be 3.1. 

We also estimated black hole masses for the AGN in 
our subset of the Palomar sa mple using the central ve- 
locity dispersions measured by IHo et all (|2009h . We con- 
verted velocity dispersion to black hole mass using the 
relation: logfMBH/Mfl ) = 8.13 + 4.02 log(cr/200kms"^) 
iTremaine et al.1 120021 '). For consistency this method was 
used for all of the galaxies, including those for which direct 
kinematic measurements were also available. 



2.2 HST data reduction 

HST images of the galaxies in our sample wer e obtained 
from the public archive. Previous studies (e.g. iRest et al] 
I2OOII ) have demonstrated that many early-type gala^xies con- 
tain dust structures that may affect the brightness profiles. 
It is especially important to deal with the effects of dust in 
the AGN sample we study he re, since dust structures are al - 
most ubiquitous in LLAGN (jConzalez Delgado et al.ll2008l ) 



and are probably more comm on in LLAGN than in normal 
galaxies l|Martini et al.|[2003l ). Therefore to minimise the ef- 
fects of dust we preferred to use images taken in the near- 
infrared using NICMOS, with the F160W filter. However, for 
many of the galaxies in our sample NICMOS images were 
unavailable, so for these we used images taken in the optical 
using WFPC2 or ACS. Only broadband filters were used, 
and preference was given to filters at longer wavelengths — 
any shorter than F555W were rejected. Using different fil- 
ters might affect our results if the brightness profile depends 
on the wavelengths at whic h the galaxy is observed, however 
iRavindranath et al.l (|200ll ) showed that the classification of 
a galaxy as core or power-law, based on Nuker fits, is in- 
dependent of the wavelength that is used. We were able to 
obtain suitable images for 110 galaxies in our sample; the 
basic properties of these galaxies are presented in table [1] 
As described in Section 12.11 the data p resented in this ta- 
ble were obtained from Ho e t al] (| 19971 ). although data was 
only available for 107 of the 110 Palomar survey galaxies for 
which we were able to obtain HST images. 

The images from WFPC2 and ACS were calibrated us- 
ing the standard On The Fly Reprocessing (OTFR) sys- 
tenSfl however the NICMOS images calibrated in this way 
contained a "pedestal" , which is a time- variable bias that is 
different in each quadrant. It is believed that this is caused 
by small temperature changes, which are different in each 
quadrant because they use separate amplifiers. To correct 
for the pedestal the standard CalnicA software that is used 
in the NICMOS OTFR system was applied using IRAF to 
the raw images, omitting all of the steps after and including 
flat field correction. The task BIASEQ in IRAF was then 
used to equalise the different biases in the 4 quadrants. Fi- 
nally, CalnicA was again applied using all of the steps after 
and including flat field correction, and omitting the steps 
that had already been performed. The final images still con- 
tain a bias, but this is now constant over the whole image 
and so can be accounted for by including a constant back- 
ground level when fitting the profile model to the image. 
These NICMOS images also contain several bad pixels — 
some of these are cold pixels, which appear dark because 
they have an abnormally low sensitivity; and some of these 
are hot pixels, which appear bright because they have an 
abnormally high dark current. There are also regions, called 
grots, that appear dark because flecks of antireflective paint 
scraped off of the optical baffles have reduced their sensi- 
tivity. These pixels are identified in the data quality files 
and are masked in the following analysis. The NIC2 cam- 
era also has a coronographic hole for observing targets that 
are close to bright objects, but when it is not used it still 
appears in the image as a bright patch, caused by emission 
from warmer structures that are behind it. This hole moves 
over time, which results in 2 patches in the image due to us- 
ing reference files that were taken at a different time. These 
regions have been masked. 



^ http://www.stsci.edu/instruments/wfpc2/Wfpc2_ 
dlib/WFPC2Jongdhbcover.html 

^ http:/ /www. stsci.edu/hst/acs/documcnts/handbooks/ 
DataHandbookvS/ ACS_longdhbcover.html 
^ http:/ /www. stsci.edu/hst/nicmos/documents/handbooks/ 
DataHandbookv7 / 
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Table 1. Basic properties of the galaxies in our sample for which HST images were obtained. 
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From left to right, columns show: Galaxy name; numerical Hubble type index; distance (Mpc); total B-band absolute magnitude 
of galaxy; bulge B-band absolute magnitude of galaxy; log black hole mass; log [Olll] luminosity (erg s~^); log H/3 luminosity 
(erg s~^); log 15 GHz radio luminosity (erg s~^); AGN optical spectral classification (S:Seyfert, LiLINER, Titransition, numbers 
denote type 1/2/i ntermediate, m ultiple classifications denote intermediate types, : and denote uncertain and highly uncertain 
classifications. See lHo et al.l l ll997i) for details); nuclear surface brightness profile classification, based on Sersic, Core-Sersic and 
Double-Sersic fits (see text). For brightness profile classifications with a question mark, the fits were uncertain. 
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Table 1 - continued Basic properties of the galaxies in our sample for which HST images were obtained. 
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3 PROFILE FITTING 



3.1 Galaxy Models 



We fitted the brightne ss profi l es of a ll galaxies in our sample 
with two models. The ISersi3 (|l968l ) model is parameterised 



as follows: 



I{r) = le exp -bn 



l/n 



(1) 



This model is described by three parameters: the effective 
radius r^, the intensity at the effective radius le, and the 
Sersic index n, which parameterises the curvature of the 
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profile. For 1 ^ n ^ 10, the quantity bn can be approximated 
by 6„ ^ 1.9992n - 0.3271 iCaon et al.lll993l). We a lso us e 
the Core-Sersic model introduced bv lcraham et al.l l|2003h . 
which combines a Sersic profile at large radii with a power 
law near the centre: 



an 7/0 



exp 



l/ari 



(2) 



This model uses six parameters, including the Sersic index 
n and effective radius Ve of the outer Sersic profile. The 
inner power law profile has a logarithmic slope 7, and the 
break radius rj, is the radius separating the inner power law 
and outer Sersic profiles. The parameter a determines the 
sharpness of the transition between the two profiles, with a 
high value of a indicating a sharp transition. The intensity /' 
is related to the intensity at the break radius, /;,, as follows: 



exp 



l/n' 



(3) 



The quantity 6„ is calculated as in the Sersic model. 

As noted in section 2.1 we excluded late type spiral 
galaxies from our sample as their brightness profiles are 
dominated by the disc, however there are still several early 
type spiral galaxies in our sample. It is possible that there 
is still a significant contribution from the disc component in 
the central surface brightness profiles of these galaxies, so 
for all of the spiral and lenticular galaxies in our sample we 
also added an exponential disc component to the model. 

The profiles of several galaxies in our sample showed 
an unresolved nuclear source due to the AGN. For these 
galaxies a central point source was added to the Sersic and 
Core-Sersic models. 

C ote et al.. (.2007h demonstrate that low and intermedi- 
ate luminosity early-type galaxies often contain a resolved 
nuclear source, for example a compact spheroidal or flat- 
tened stellar component. They model the brightness profiles 
of these galaxies using a Double-Sersic model: 



/(r) =/e,i exp 



-I- Ie,2 exp 



1/nl 



l/n2 



(4) 



The inner Sersic model (parameterised by 7e,i, re,i and nl) 
describes the resolved nuclear component while the outer 
Sersic profile (parameterised by 7e,2, re,2 and n2) describes 
the host galaxy. The brightness profiles of several galaxies 
in our sample showed a central excess above a single Sersic 
profile that was too extended to be described by a point 
source, so these galaxies were also fitted with this Double 
Sersic model. 

Using the fits with the Sersic and Core-Sersic laws, we 
classified a galaxy as 'Core' i f it satisfied the fol lowing cri- 
teria (based on those used bv lTruiillo et al.lbOO'jl '): 

(i) The Core-Sersic fit mus t give a lower r educe d than 
the Sersic fit. We note that iTruiiUo et al] (|2004l ') required 
that the reduced of the Core-Sersic fit must be at least 
a factor of 2 less than that of the Sersic fit, however we 
found that very few galaxies in our sample satisfied such a 
strict criterion. This discrepancy is likely due to differences 
in how the data points are weighted in the fits, as they use 



equal weights whereas we weight the data points using their 
Poissonian errors, which are larger at smaller radii. 

(ii) To ensure that the core is well sampled by the data 
the break radius must be greater than the radius of the 
second data point. 

(iii) The slope of the inner power law profile, 7, must be 
less than the logarithmic slope of the best-flt Sersic model 
within the break radius. 

Galaxies for which the Double-Sersic model gave the 
lowest reduced were classified as 'Double' galaxies (i.e. 
contain a resolved nuclear source). All other galaxies were 
classified as 'Sersic' galaxies. 



3.2 Extracting ID Profiles 

The ID semi-major axis brightness profiles of each galaxy 
were extracted using the ELLIPSE task in IRAF, which fits 
elliptical isophotes to the galaxy images. Regions of dust 
that were visible in the images were masked from these fits. 
The position angle and ellipticity were free to vary with each 
isophote, enabling radial variations in these quantities, such 
as isophotal twists, to be accounted for. 



3.3 Fitting Procedure 

The Sersic, Core-Sersic and (where necessary) Double-Sersic 
models were fitted to the ID brightness profiles by minimis- 
ing the statistic using the 5/ierpt0 modelling and fitting 
software package, which is a part of the Chandra Interactive 
Analysis of Observations (CIAC0) software. The data points 
were weighted in these fits using their Poissonian errors. 

To accurately reproduce the brightness profile at small 
radii it is necessary to account for the point spread function 
(psf) of the instrument used. For each image we created a 
2D model of the psf at the position of thegalaxy's centre on 
the detector using the Tiny Tim softwarqj. It was necessary 
to create a separate psf model for each image because the 
psf depends on the position on the detector. The ELLIPSE 
task in IRAF was used to extract the ID profile of the psf, 
then the galaxy model was convolved with this psf in Sherpa 
before fitting to the observed brightness profile. 

In the Core-Sersi c model the param eters a, n and 7 
tend to be degenerate. ITruiillo et alTl|2004 ) recommend that 
a should be fixed at a large value to force an abrupt transi- 
tion between the outer Sersic and inner power-law profiles, 
however we found that in several galaxies the profile was 
poorly fitted by such a fixed-a model, but it was well fitted 
if a was free to vary. Therefore we fixed q at a = 50 to 
obtain an initial fit and then we allowed a to vary freely to 
obtain the final fit. We investigated how robust the parame- 
ters were for the free-a fits and we found that the parameters 
of the outer Sersic component were uncertain in some of the 
galaxies, particularly those with a low value of a. This may 
also be caused by the limited radial extent of the fitting 
regions, which were comparable to the effective radius re 
in some galaxies. However, the uncertainties in 7 were rela- 
tively small, and the break radii rb of the Core-Sersic models 



* http://cxc.harvard.edu/sherpa/ 
^ http://cxc.harvard.edu/ciao/ 

^ http : / / www. stsci.edu /software /tinytim /tinytim. html 
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were well within the fitting region, so the classifications as 
Core or Sersic were robust. 

3.4 Results 

We were able to obtain fits for all 110 galaxies in our sample, 
however we are only confident in 62 of these fits, which are 
reported in table[21 44 of the remaining galaxies showed large 
regions of dust that made the fits uncertain, even after we 
had masked most of this dust. The profiles of the remaining 
4 galaxies could not be well described by a Sersic, Core- 
Sersic or Double-Sersic model. We report these fits in tabled 
however they will be discarded for most of the following 
analysis. We present examples of confident fits using the 
Sersic, Core-Sersic and Double-Sersic models in fig. [1] and 
two uncertain fits in fig. [2] showing for each the galaxy image 
and the radial profiles of the galaxy and the best fitting 
model. 

We included an exponential disc component in the mod- 
els for the spiral and lenticular galaxies in our sample, how- 
ever in many cases this component was negligible. We there- 
fore only include the disc component parameters in tables [2] 
and O for the galaxies where it is significant. 

To enab le us to co mpare our results with previous stud- 
ies such as iCapetti fc B almaverdc (2005), which used the 
Nuker classification scheme, we also fitted the Nuker model 
to all galaxies in our sample. If we assume that 'intermedi- 
ate' as well as 'core' galaxies from the Nuker scheme cor- 
respond to 'Core-Sersic' galaxies and 'Double-Sersic' and 
'Sersic' galaxies both correspond to 'power-law' galaxies in 
the Nuker scheme then the Nuker classifications of the galax- 
ies in our sample agree with those that we present in table [5] 
for 82% of the galaxies for which we have a confident fit. 
This suggests that while the classification of many galaxies 
in our sample remain unchanged when we use the Nuker 
model, there are still some discrepancies between the two 
classification schemes. 



4 DISCUSSION 

Of the 62 galaxies for which we are confident in the bright- 
ness profile fits, we have 21 Core galaxies, 26 Sersic galax- 
ies and 15 Double-Sersic galaxies. We only have both lu- 
minosity data and a confident classification for 60 of the 
197 galaxies in the Palomar sample that are confirmed as 
AGN hosts, so we need to verify that we have not intro- 
duced any observational biases, for example if bright galax- 
ies are more likely to be imaged with //ST then there would 
be a disproportionately large number of them in our final 
sample. In fig. Owe compare the distributions of the Hub- 
ble type, bulge B band magnitude, nuclear radio luminos- 
ity and [Olll] line luminosity for the 60 galaxies in our fi- 
nal sample, the 107 galaxies for which we have a confident 
or uncertain fit and luminosity data, and the 197 galax- 
ies ho sting AGN in the Palomar sample, that were stud- 
ied in lNagar et alj l|2005l ) (the Nagar sample). Wc used the 
ASURV Rev 1.2 software (|LaVallev. Isobe fc Fei gclson 1992) 
to test the significance of correlations and differences be- 
tween distributions for censored data, which includes non- 
detections with li mits. This software i mplem ents the meth- 
ods presented in iFeigelson" fc NelsonI (| 19851 ) for univariate 



problems and llsobe. Feigelson fc NelsonI l| 19861 ) for bivariate 
problems. 

According to the generalised Wilcoxon tests the distri- 
butions of Hubble type in our final sample and the Nagar 
sample are inconsistent with being drawn from the same dis- 
tribution at the 99.99 per cent confidence level. We can see in 
fig. [S] top left panel, that our final sample has a higher pro- 
portion of elliptical galaxies, which is to be expected because 
we discarded all late spirals, and a disproportionately high 
number of the fits for early spirals were uncertain because 
they are more likely to contain dust. Using the generalised 
Wilcoxon tests, we also find that the distributions of bulge 
magnitude in our final sample and the Nagar sample are only 
marginally inconsistent with being drawn from the same dis- 
tribution at the 90 per cent confidence level. There are a 
small number of galaxies in the Nagar sample with very dim 
bulges that we missed because we discarded late type spirals. 
Also, we were more likely to find HST images for brighter 
galaxies. Of the galaxies for which we did obtain HST data, 
a high proportion of those with intermediate bulge magni- 
tudes had uncertain Nuker fits, probably because these tend 
to correspond to later type galaxies which are more likely to 
be dusty. According to the generalised Wilcoxon tests the 
distributions of nuclear radio luminosities in our final sam- 
ple and the Nagar sample are also marginally inconsistent 
with being drawn from the same distribution at the 90 per 
cent confidence level. From fig. O bottom left panel, we can 
see that there is a higher proportion of the small number 
of AGNs with higher radio luminosities in our final sample. 
Finally, using the generalised Wilcoxon tests we find that 
the distributions of [Olll] line luminosities in our final sam- 
ple and the Nagar sample are consistent with being drawn 
from the same distribution, although we can see in fig. O 
bottom right panel, that the galaxies with the lowest [Olll] 
line luminosities do not appear in our final sample. 

4.1 Double-Sersic Galaxies 

We find 15 galaxies that are best fit by a Double-Sersic 
model. Previous studies (e.g. ICote et"al] [2007 ) attributed 
the inner Sersic component in such models to a resolved 
stellar nucleus, with a radius « 0.02re, where re is the ef- 
fective radius of the outer component. However, several of 
the Double-Sersic galaxies in our sample have inner compo- 
nents with a much larger radius. This suggests that they 
might not be due to resolved stellar nuclei, for example 
some of them may ins tead be inner discs or nuclear bars 
(|Erwin fc Sparkel [20021 ). We studied the isophotes of the 
galaxies and we did not find any with flattened isophotes 
towards the centre, which suggests that nuclear bars are 
unlikely. However, searching through the literature we flnd 
that inner discs have been discovered in som e of the galax- 
ies in our sample, for example NGC 3945 ([Moiseev et al.[ 
|2004 ) and NGC 7217 l|Zasov fc Sirchenkd ll997t). We there- 
fore cannot be certain whether the Double-Sersic galaxies in 
our sample contain resolved stellar nuclei, or if they instead 
contain other nuclear structures such as inner discs. 

4.2 Black hole mass and bulge magnitude 

Fig. [4l left panel, shows the distributions of central black 
hole mass for Core, Sersic and Double-Sersic galaxies. The 
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Figure 2. Galaxy images (top row) and radial brightness profiles (bottom row) for a fit that is uncertain due to dust (IC 356; left) and 
a fit that is uncertain because the galaxy has a complex profile (NGC 4111; right). 
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Table 2. Best fit parameters of the Core, Sersic and Double Sersic galaxies from the confident fits. 
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6.06 


0.74 


0.16 


8.86 


215.30 


3.88 


- 


NGC7626 


NICMOS 


F160W 


9.26 


0.47 


0.30 


5.39 


19.73 


2.79 


- 








Sersic Galaxies 










Galaxy 


Instrument 


Filter 




n 


rl 


Me, disc 












F 1 ROW 

r _LUU VV 


18.12 


2.09 


4.22 






90 'i7 






NICMOS 


FT fiflW 

C _LUU VV 


20.26 


3.53 


22.83 






22.05 






WFPP9 
VV r 1 v_/ z 


FS1 AWJ 
r o i-'-t VV 


20.93 


3.91 


29.76 














FT ROW 

C ±UU VV 


19.50 


3.45 


16.52 






20.77 








FT ROW 
r J.UU VV 


17.80 


4.07 


3.63 






1 7 iiQ 




NGG'^94^ 


i N ± v_y iV J. V_/ O 


FT ROW 

C -LUU VV 


17.59 


3.03 


6.02 






19.52 






VV r 1 v_,'Z 


FS1 AW 

P 1^ VV 


22.38 


5.95 


48.81 














FT ROW 

-T _LUU VV 


16.64 


1.43 


2.24 






17.17 








F 1 ROW 

r _LUU VV 


15.95 


1.50 


1.04 


-LO.UO 


A 


19 21 




IN vjv^oyuu 


NTPMOS 

i > ± \_' iV J. V, / o 


Fl ROW 
S: _LUU VV 


22.22 


6.55 


38.11 














F 1 ROW 

r _LUU VV 


19.60 


2.47 


4.84 






9n 1 7 

ZU. J- 1 








FT ROW 

C _LUU VV 


18.21 


3.10 


4.52 














FT ROW 

r _LUU VV 


18.76 


4.30 


9.00 














r iuU VV 


16.76 


1.93 


5.27 




ooy. / o 


1 CI 




NGC4203 


ACS 


F814W 


21.49 


4.30 


22.92 






21.22 




NGC4314 


NICMOS 


F160W 


22.90 


5.41 


154.78 










NGC4378 


WFPC2 


F606W 


19.75 


2.78 


5.93 


21.51 


24.46 






NGC4450 


WFPC2 


F814W 


25.61 


8.01 


502.69 










NGC4459 


NICMOS 


F160W 


18.12 


3.56 


8.62 










NGC4636 


NICMOS 


F160W 


17.87 


1.20 


7.21 










NGC4698 


WFPC2 


F814W 


16.65 


2.03 


1.64 


18.81 


11.65 


20.60 




NGC5273 


NICMOS 


F160W 


17.85 


1.66 


2.08 






19.90 




NGC5838 


NICMOS 


F160W 


20.50 


5.21 


37.30 










NGC5982 


NICMOS 


F160W 


17.86 


2.03 


5.54 










NGC7742 


NICMOS 


F160W 


19.94 


4.64 


9.55 










NGC7743 


NICMOS 


F160W 


20.17 


6.34 


9.66 











Corrected for extinction using values from IHo et ah and converted to the Johnson V band, assuming that the 

galaxy has the spectrum of a K2 giant. 
Measured in arcsec. 

Point source magnitude, corrected for extinction using values from lHo et al. I lll997h and converted to the Johnson V band, 
assuming a spectral index of 1. 
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Table 2 - continued Best fit parameters of the Core, Sersic and Double Sersic galaxies from the confident fits. 



Double Sersic Galaxies 



Galaxy 


Instrument 


Filter 


Mel 


ni 


^el 


Me2 


n2 


^e2 








NGC404 


NICMOS 


F160W 


19.74 


15.72 


2.55 


24.97 


6 


35 


712.26 


~ 


~ 




NGC3169 


NICMOS 


F160W 


14.90 


0.81 


0.38 


17.24 


1 


60 


6.72 








NGC3945 


WFPC2 


F814W 


17.12 


1.57 


1.10 


18.67 


1 


04 


9.33 






20.60 


NGC3998 


ACS 


F814W 


17.11 


3.04 


0.35 


19.66 


2 


90 


9.75 








NGC4125 


NICMOS 


F160W 


16.55 


1.04 


0.84 


18.46 


2 


32 


17.56 








NGC4143 


NICMOS 


F160W 


18.24 


0.48 


3.25 


17.69 


2 


51 


5.59 






21.12 


NGC4435 


NICMOS 


F160W 


16.12 


0.59 


0.39 


16.94 


1 


64 


5.52 






19.56 


NGC4550 


NICMOS 


F160W 


16.96 


0.16 


1.92 


19.12 


3 


56 


5.10 








NGC4725 


NICMOS 


F160W 


15.20 


1.06 


0.13 


18.70 


2 


83 


14.10 








NGC4736 


NICMOS 


F160W 


16.14 


0.56 


8.25 


17.98 


3 


54 


19.42 








NGC5377 


NICMOS 


F160W 


15.06 


0.95 


0.26 


16.91 





63 


2.81 


20.41 


440.82 




NGC5678 


NICMOS 


F160W 


15.68 


0.16 


0.02 


23.47 


6 


07 


133.73 


22.11 


367.00 




NGC5746 


NICMOS 


F160W 


17.44 


0.16 


1.41 


17.56 


2 


09 


6.74 


19.06 


623.32 


20.61 


NGC6703 


WFPC2 


F814W 


20.57 


0.82 


10.87 


21.95 


6 


22 


25.78 








NGC7217 


NICMOS 


F160W 


16.52 


0.93 


0.79 


19.60 


2 


64 


20.05 


20.51 


167.78 





5 
Hubble Type 



LogL^_, 
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-15 



-20 
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Figure 3. Histograms showing the distributions of Hubble type (top left), bulge B band magnitude (top right), nuclear radio luminosity 
(bottom left) and [Oiii] line luminosity (bottom right) in our final sample (dotted line), the galaxies for which we have a confident or 
uncertain fit and luminosity data (dashed line) and the Nagar sample (solid line). The arrows denote numbers of galaxies in the histogram 
bins which are limits only. 
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Table 3. Best fit parameters of the Core, Sersic and Double Sersic galaxies from the uncertain fits. 













Core Galaxies 














Galaxy 


Instrument 


Filter 






7 


n 


j.b 
' e 


a 




Uncertainty 






i > VJ V_y> J _L \J 


WFPC2 


F814W 


12.51 


0.89 


0.00 


5.06 


100.29 


3.27 


21.27 


D 






VJ V_y 'J (J'-t\J 


WFPC2 


F702W 


11.23 


1.24 


0.00 


6.54 


419.02 


4.40 




D^'^t^ 
us y 
















Sersic Galaxies 














Galaxy 


Instrument 


Filter 


MS 


n 




f^e,disc 






Uncertainty 








IC356 


WFPC2 


F814W 


25.83 


9.16 


1491.05 


21.32 


2904.34 


- 


Dusty 








NGC660 


ACS 


F814W 


19.25 


0.23 


8.98 


_ 


- 


- 


Dusty 








NGC2273 


NICMOS 


F160W 


16.17 


0.74 


1.82 


19.55 


15.45 


18.76 


Dusty 








NGC2681 


NICMOS 


F160W 


23.12 


12.35 


103.2 


_ 


- 


- 


Complex 








NGC2683 


NICMOS 


F160W 


22.79 


6.17 


184.16 


19.43 


369.32 


- 


Dusty 








NGC3190 


ACS 


F814W 


18.85 


1.97 


9.54 


_ 


- 


- 


Dusty 








NGC3368 


NICMOS 


F160W 


17.89 


2.42 


14.26 


_ 


- 


16.81 


Dusty 








NGC3489 


WFPC2 


F814W 


19.51 


4.22 


15.24 


_ 


- 


- 


Dusty 








NGC3507 


WFPC2 


F606W 


21.01 


4.8 


7.95 


21.75 


43.06 


- 


Dusty 








NGC3623 


WFPC2 


F814W 


21.16 


5.81 


63.13 


19.68 


579.52 


- 


Dusty 








NGC3626 


NICMOS 


F160W 


16.67 


2.32 


3.14 


_ 


- 


17.87 


Dusty 








NGC3627 


NICMOS 


F160W 


17.79 


2.78 


10.46 


19.57 


344.24 


- 


Dusty 








NGC3628 


NICMOS 


F160W 


25.21 


5.6 


5850.16 


_ 


- 


- 


Dusty 








NGC4013 


WFPC2 


F814W 


19.37 


0.18 


7.9 


20.55 


228.1 


- 


Dusty 








NGC4281 


WFPC2 


F606W 


25.97 


8.71 


509.97 




- 


- 


Dusty 








NGC4388 


NICMOS 


F160W 


24.09 


7.79 


121.02 




- 


- 


Dusty 








NGC4429 


WFPC2 


F606W 


22.78 


4.52 


134.39 




- 


17.43 


Dusty 








NGC4438 


WFPC2 


F814W 


18.58 


3.42 


12.99 




- 


- 


Dusty 








NGC4477 


WFPC2 


F606W 


21.3 


3.58 


35.93 




- 


- 


Dusty 








NGC4494 


WFPC2 


F814W 


33.22 


17.85 


46101.91 








Dusty 








NGC4565 


NICMOS 


F160W 


16 


1.44 


3.9 


18.3 


61.7 


18.66 


Dusty 








NGC4579 


WFPC2 


F791W 


19.24 


3.32 


9.77 


20.22 


48.67 


19.39 


Dusty 








NGC4772 


WFPC2 


F606W 


23.7 


6.27 


68.5 


21.68 


36.77 


- 


Dusty 








NGC4826 


NICMOS 


F160W 


16.68 


2.21 


6.87 


19.5 


80.47 




Dusty 








NGC4S66 


ACS 


F814W 


33.20 


13.65 


33227.88 








Dusty 








NGC,'S44S 


NICMOS 


J. _L \J\J V V 


17.36 


1.49 


2.56 


20.71 


130.06 












NGCB.'ififi 


WFPC2 


FfiOfiW 


17.16 


0.86 


5.23 








D 








NGCSSfifi 


NICMOS 


F1 fiflW 


18.26 


1.77 


23.21 








Dusty 








NGC,5q8,'i 

1 ^ V_J V_y coo 


NICMOS 


FlBtlW 

J. _L \J\J V V 


19.56 


2.03 


4.91 






22.46 


Complex 








NGCfi34n 


ACS 


F814W 


20.28 


3.44 


5.07 


22.14 


27.77 




Dusty 








NGC7177 


NICMOS 


J. J_ \J\J V V 


18.16 


1.89 


5.97 
















NGC7331 


NICMOS 


Fl fiflW 


22.96 


8.09 


344.32 


18 


21.64 












NG07814 


NICMOS 


J. _L \J\J V V 


17.82 


1.91 


13.16 








Dusty 
















Double Sersic Galaxies 














Galaxy 


Instrument 


Filter 




ni 


' el 






' e2 








Uncertainty 


NGC2768 


WFPC2 


F814W 


17.14 


1.00 


0.44 


24.88 


5.98 


547.44 








Dusty 


NGC2787 


WFPC2 


F814W 


16.74 


1.54 


0.58 


19.78 


2.20 


20.34 








Dusty 


NGC2859 


ACS 


F814W 


20.00 


0.16 


0.31 


20.20 


2.32 


9.46 








Complex 


NGC3675 


NICMOS 


F160W 


19.14 


0.17 


8.96 


21.11 


5.51 


42.71 


21.11 


1141.66 




Dusty 


NGC3705 


WFPC2 


F814W 


14.32 


0.4 


0.22 


20.58 


4.01 


28.64 


20.82 


215.19 




Dusty 


NGC4111 


NICMOS 


F160W 


17.58 


0.16 


1.41 


18.75 


3.69 


19.13 








Complex 


NGC4220 


ACS 


F814W 


20.60 


2.25 


5.57 


20.29 


0.14 


46.12 








Dusty 


NGC4293 


NICMOS 


F160W 


16.4 


0.43 


0.82 


29.77 


8.65 


25512.11 


22.23 


282.46 


17.99 


Dusty 


NGC4394 


ACS 


F814W 


19.83 


0.16 


1.02 


21.45 


3.81 


9.06 


22.03 


24.62 




Dusty 


NGC4501 


WFPC2 


F606W 


15.97 


0.73 


0.64 


19.06 


1.11 


12.09 


22.08 


202.65 


19.89 


Dusty 


NGC4569 


NICMOS 


F160W 


15.15 


3.72 


0.72 


15.86 


0.59 


2.36 


18.95 


22.24 


18 


Dusty 


NGC4596 


WFPC2 


F606W 


15.54 


0.89 


0.07 


22.45 


4.22 


66.19 








Dusty 


NGC4750 


NICMOS 


F160W 


14.46 


0.16 


0.07 


17.6 


1.36 


2.39 


20.36 


157.09 




Dusty 



Corrected for extinction using values from iHo et ah and converted to the Johnson V band, assuming that the 

galaxy has the spectrum of a K2 giant. 
^ Measured in arcsec. 

Point source magnitude, corrected for extinction using values from I Ho et al.l l|l997tl and converted to the Johnson V band, 
assuming a spectral index of 1. 
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6 7 8 9 

Log M3„/M^^, 

Figure 4. Histograms showing the distributions of central black hole 
(top), Sersic galaxies (middle) and Double-Sersic galaxies (bottom). 



black hole masses of the Core galaxies are typically higher 
than those of the Sersic or Double-Sersic galaxies, for exam- 
ple the median of log Mbh/Mq is 8.39 for the Core galaxies, 
compared to 7.57 and 7.84 for the Sersic and Double-Sersic 
galaxies respectively. We also find that the most massive 
black holes are almost exclusively found in Core galaxies — 
with the exception of one Double-Sersic galaxy, NGC 3998, 
all SMBHs with log Mbh/Mq > 8.5 are found in Core galax- 
ies. Conversely, all SMBHs with log Mbh/Mq < 7.6 are 
found in Sersic or Double- Sersic galaxies. This is consis tent 
with previous studies (e.g. ICapetti fc Balmaverd3l2006l '). 

Fig. 21 right panel, shows the distributions of bulge B 
band magnitude for Core, Sersic and Double-Sersic galaxies. 
The bulges of Core galaxies are generally brighter than those 
of Sersic or Double-Sersic galaxies, although they overlap 
significantly — for all Core galaxies AlB,buige ^ —18.6, and 
for all Sersic galaxies MB^buigc ^ —20.4. Most Double-Sersic 
galaxies are within the same range as the Sersic galaxies, 
although there are two examples t hat are more luminous 
and o ne that is much less luminous. ICapetti fc Balmaverdd 
found a similar trend in K band magnitudes, although 
they also found core galaxies covering a much broader range 
of magnitudes. It is also worth noting that, with the ex- 
ception of one Double-Sersic galaxy, NGC 4125, all galaxies 
with bulge B magnitudes below -21 are Core galaxies. 

We note that the most massive black holes being found 
in Core galaxies is to be expected as they are at the more 
lumino us end of the luminosity function (e.g. lFerrarese et al.l 
l2006bl . and also confirmed for our sample by Figure [H right 
panel), so the y also tend to have larger velocity dispersions 
IjFaber fc Jac kson 1976) and hence black hole masses, which 
we derived from the velocity dispersions. 

4.3 Radio-Loudness 

The radio-loudness of an AGN is typically measured as the 
ratio of the radio to optical luminosity of the nucleus, how- 
ever we were only able to extract a point source in 23 of the 
galaxies for which we are confident in the Core/Sersic fits, 
so this would limit the number of nuclei for which we can 
measure the radio-loudness. Furthermore this could bias any 
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Figure 5. Histograms showing the distribution of iiy,ra[i/^[Olll] 
for Core galaxies (top), Sersic galaxies (middle) and Double-Sersic 
galaxies (bottom) . The arrows denote numbers of galaxies in the 
histogram bins which are limits only. 

correlations that we observe because we would be preferen- 
tially selecting bright, type 1 nuclei that are unobscured by 
nuclear dust. We therefore use the narrow [Olll] line lumi- 
nosity instead of the total optical luminosity of the nucleus, 
because we have this data for most of the galaxies in our 
sample, and since [Olll] line emission originates from the 
narrow line region, there should not be a bias between type 
1 and type 2 nuclei. 

In Fig. [5] we show the distribution of the radio-loudness 
[Oiii] for Core, Sersic and Double-Sersic galaxies. 
According to the generalised Wilcoxon tests the distribu- 
tions of radio-loudness in Core and Sersic galaxies are in- 
consistent with being drawn from the same distribution at 
the 99.9 per cent confidence level, while the distributions of 
radio-loudness in Sersic and Double-Sersic galaxies are con- 
sistent with being drawn from the same distribution. This 
confirms that the radio-loudness of an AGN is related to 
whether its host galaxy is classified as a Core galaxy. We 
can see in fig. [5] that the most radio-loud AGNs are hosted 
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in Core galaxies, however AGNs with \og{L^^radio/ Lyoiii]) ^ 
—0.8 can be found in both Core and Sersic galaxies. This 
large ove rlap has not been obser v ed in previous studies, for 
example ICapetti fc Balmaverdd (|2006l ) found that almost 
all core galaxies were more radio-loud than the power-law 
galaxies (classified using the Nuker scheme) , with very little 
overlap. Their sample was restricted to early-type galaxies 
whereas our sample also included early spiral galaxies, how- 
ever excluding the spirals in our sample did not recover the 
obvious split that they observed. The optical nu clear lumi- 
nosities reported in lCapetti fc Balmaverdj (|2006l ) are higher 
for power-law galaxies than for core galaxies, however the 
lowest [Olll] line luminosities of power-law galaxies are as 
low as those of core galaxies, although they do extend to 
higher [Olll] line luminosities than core galaxies. This sug- 
gests that the discrepancy between the correlation that we 
observe and that observed by Capetti & Balmavcrde (200^) 
is due to the different definitions used for the radio-loudness. 
To test this possibility we looked at the galaxies for which we 
were able to extract a point source in the profile fitting and 
we used this point source to determine the ratio of nuclear 
radio luminosity to nuclear optical continuum luminosity, 
Lv,radio/Lo. FoT these galaxies we then compared the dis- 
tributions of both L^^radio/Lo and L^^radio/ L\^oiii^ for Core, 
Sersic and Double-Sersic galaxies; these results are presented 
in Fig. [6] According to the generalised Wilcoxon tests the 
distributions of L^j^radio/ Lo in Core and Sersic galaxies are 
inconsistent with being drawn from the same distribution at 
the 99.99 per cent confidence level, while the distributions 
of Lij^radio/ Lyoiii] ^rc incousisteut at the 99.9 per cent level. 
Furthermore, we can see from Fig.[6]that there is no overlap 
between the Core and Sersic subsamples when we use the nu- 
clear optical continuum luminosity. This demonstrates that 
the split between Core and Sersic galaxies is clearer when 

we use Lv,radiolLo- 

We also looked at using the narrow H/3 line luminos- 
ity instead of [Olll] . Fig. [7] shows that using this definition 
for the radio-loudness produces the same difference between 
Core and Sersic galaxy radio-loudness distributions that was 
observed when [Olll] was used. 

4.4 Correlations of radio-loudness with black hole 
mass and bulge magnitude 

We saw in § 14.21 that the distributions of both the black 
hole mass and the bulge B band magnitude are different for 
Core and Sersic galaxies, so this could lead to the observed 
trends of radio-loudness with the Core/Sersic classification 
if radio-loudness is dependent on one of these properties. To 
test this possibility we plot radio-loudness against the black 
hole mass and the bulge B band magnitude; these plots are 
presented in Fig. [S] 

Using Spearman's rho to test the correlation between 
radio-loudness and black hole mass, we find that there is 
a correlation at the 99.99 per cent confidence level, with 
galaxies hosting high mass black holes generally being more 
radio-loud. Similar correlations have been demonstrated by 
previous studies of AGN, for example radio-loud quasars 
are found to contain system atically more massive black 
holes than radio-quiet quasa rs l|Laoil(2000l : lMcLure fc Jarvid 
|2004 ). and iDoi et"al] l|2006l ) found a best-fit plane in the 
three dimensional space of radio-loudness, black hole mass 




Figure 7. Histograms showing the distribution of Lv^radl^np for 
Core (top), Sersic (middle) and Double-Sersic (bottom) galaxies. 
The arrows denote numbers of galaxies in the histogram bins 
which are limits only. 



and Eddington ratio for a sample of 48 LLAGN taken from 
the Palomar spectroscopic survey. In fig. [S] left panel, we 
can see that there is still considerable scatter between these 
variables, which suggests that other factors may be involved. 

We also used Spearman's rho to test the correlation be- 
tween radio-loudness and bulge B band magnitude, and we 
found that there is a correlation at the 99 per cent con- 
fidence level, with luminous galaxies generally being more 
radio loud. However, fig. [H right panel shows that there is 
also scatter between these two variables, which again sug- 
gests that there may be other factors involved. 

It is possible that the scatter in the correlations of radio- 
loudness with bulge magnitude and black hole mass can be 
explained by the Core/Sersic/Double-Sersic classification. 
To test correlations with this classification we parameterised 
the galaxies using the central deficit or excess in the surface 
brightness profile compared to the inward extrapolation of 
the outer S ersic component, using the A0.02 parameter in- 
troduced bv lCote et all (|2007l '). This parameter is defined as 
A0.02 = log Cg / Cs , where Cg is the total luminosity of the 
best fit galaxy model within a radius of 0.02re and Cs is 
the total luminosity of the outer Sersic component in this 
region. Core galaxies have A0.02 < 0, Sersic galaxies have 
A0.02 = and Double-Sersic galaxies have A0.02 > 0. We 
then tested the partial correlations for both black hole mass 
and bulge magnitude with radio loudness and A0.02, using 
the cens_tau code H, wh ich uses the methods presented in 
lAkritas fc SiebertI ifl996t ) to search for partial correlations in 
censored data. These partial correlations test whether there 
are residual correlations between radio-loudness and A0.02 
after the overall correlations have been accounted for. 

Zero partial correlation between the bulge magnitude, 
radio-loudness and A0.02 is rejected at the 99 per cent con- 
fidence level, and zero partial correlation between the black 
hole mass, radio-loudness and A0.02 is rejected at the 99.99 
per cent confidence level. These results suggest that there is 
a significant residual correlation between the radio-loudness 



http: / /astrostatistics. psu.edu/statcodes/cens_tau 
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Figure 6. Histograms comparing the distributions of L^ j-adio/^o (left panel) and ii/,ra(iio/^[Olll] (right panel) for the Core (top), 
Sersic (middle) and Double-Sersic (bottom) galaxies for which we were able to extract a point source. The arrows denote numbers of 
galaxies in the histogram bins which arc limits only. 
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Figure 8. Radio-loudness i^,, ad/^[Olll] plotted against black hole mass (left panel) and against bulge B band magnitude 
Core galaxies are represented by filled squares, Sersic galaxies by empty circles and Double-Sersic galaxies by empty triani 



(right panel), 
gles. 



of an AGN and the classification of its host galaxy as Core, 
Sersic or Double-Sersic, after the effects of bulge magni- 
tude and black hole mass have been accounted for. How- 
ever, iDoi et al.l l|2006h demonstrated that the radio-loudness 
correlates with both the black hole mass and the accre- 
tion rate in a sample of 48 LLAGN, so we also need to 
consider whether A0.02 correlates with the accretion rate. 
We used the ratio of [Olll] line luminosity to Eddington 
luminosity, log{L[oiii]/LEdd), as a proxy for the accre- 
tion rate and tested the partial correlation of A0.02 with 
\og{L[oiii]/ Lsdd) and black hole mass, following the same 
method as above. We found that A0.02 is partially corre- 
lated with \og{L[oiii]/ Lsdd) at the 99.99 per cent confi- 
dence level. Therefore it is possible that the partial correla- 
tion between radio-loudness and A0.02 arises because Core 
galaxies tend to have a higher central black hole mass and a 
lower accretion rate. This possibility needs further inves- 
tigation. The correlation between accretion rate and the 
classification of the host galaxy as Core, Sersic or Double- 
Sersic is illustrated in fig.lH] where we plot log L^^rad/ Lyoiii] 



against logLyoiii]/ Lsdd for the Core, Sersic and Double- 
Sersic galaxies. 



5 CONCLUSIONS 

In this paper we have studied the brightness profiles of el- 
liptical and early type spiral galaxies in the Palomar spec- 
troscopic survey that are confirmed as AGN hosts. We fit- 
ted Sersic, Core-Sersic and, where necessary, Double-Sersic 
models, plus a point source where needed, to the ID semi- 
major axis brightness profiles extracted from high resolution 
HST images of these galaxies. By comparing these fits we 
were able to classify the galaxies as Core, Sersic or Double- 
Sersic galaxies, and then we investigated how other proper- 
ties of the host galaxy and the AGN relate to this classifi- 
cation. 

We found that Core galaxies were generally more lu- 
minous and hosted higher mass black holes than Sersic 
and Double-Sersic galaxies, although there was considerable 
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Figure 9. Radio-loudness L^ ^-^d/ Lyoiii] plotted against the 
[Olll] line luminosity as a fraction of the Eddington luminosity, 
^[OIII]/ ^Edd) for Core galaxies (filled squares), Sersic galaxies 
(empty circles) and Double-Sersic galaxies (empty triangles). 



overlap between thes e subsamples. These r esults agree with 
previous studies (e.g. lCapetti fc Balmaverda,200& . who used 
the Nuker classification scheme). 

To measure the radio-loudness of the AGNs we took the 
ratio of nuclear radio luminosity to [Olll] line luminosity, be- 
cause we could only estimate the optical continuum luminos- 
ity of the AGN for approximately one third of the galaxies, 
for which we could extract a point source in the profile fit- 
ting. Furthermore, using the [Olll] line luminosity prevents 
bias against faint nuclei that are more difficult to extract. 
Using this definition of the radio-loudness we found that 
Core galaxies were generally more radio-loud than Sersic 
an d Double-Sersic galaxie s , in a greement with the results 
of ICapetti fc Balmaverdj (|2P06i) for a radio-selected sam- 
ple of AGN based on the Nuker scheme. However, we found 
significan tly more overlap betwe e n the se subsamples com- 
pared to ICapetti fc Balmaverdd (|2006[ ). This difference is 
most likely because we use a different definition of the radio 
loudness, as we find that there is much less overlap when we 
use the nuclear optical continuum luminosity. Therefore we 
conclude that the radio-loudness/brightness-profile connec- 
tion uncovered for radio-selected AGN also applies to our 
optically-selected sample. 

We also looked at how the radio-loudness, defined us- 
ing the [Olll] line luminosity, correlated with the black hole 
mass and bulge B band magnitude of the host galaxies. We 
found that, while the radio -loudness does show correlations 
with both black hole mass and bulge magnitude, which could 
explain at least some of the correlation between brightness- 
profile and radio-loudness, we still found a significant partial 
correlation with the classification of the host galaxy bright- 
ness profile as Core or Sersic. However, the host galaxy clas- 
sification is also correlated with the accretion rate, so it is 
possible that the observed connection between brightness- 
profile and radio-loudness arises because Core galaxies tend 
to have a higher black hole mass and a lower accretion rate. 
This possibility requires further investigation. 
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